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An approach to the construction of the CNOT quantum logic gate for a 4-dimensional coupled- 
qubit model with weak but otherwise arbitrary coupling has been given recently [M. R. Geller et ai, 
Phys. Rev. A 81, 012320 (2010)]. How does the resulting fidelity depend on the form of qubit-qubit 
coupling? In this paper we calculate intrinsic fidelity curves (fidelity in the absence of decoherence 
versus total gate time) for a variety of qubit-qubit interactions, including the commonly occurring 
isotropic Heisenberg and XY models, as well as randomly generated ones. For interactions not too 
close to that of the Ising model, we find that the fidelity curves do not significantly depend on the 
form of the interaction, and we calculate the resulting interaction-averaged fidelity curve for the 
non-Ising-like cases and a criterion for determining its applicability. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

The operation of a quantum computer requires the 
implementation of a universal two-qubit quantum logic 
gate, such as the CNOT The problem of constructing 
a CNOT gate has been addressed from various perspec- 
tives and for different physical systems 04111 ■ In recent 
work, Geller et al. approached the problem of CNOT 
construction from a somewhat general standpoint. Start- 
ing with a Hamiltonian for a four-dimensional coupled- 
qubit model, they derived a CNOT pulse sequence as- 
suming weak coupling. In this work, we calculate the 
intrinsic fidelity of the CNOTs constructed according to 
this protocol. 

The question of CNOT gate fidelity has already been 
discussed from other standpoints [3 . Fidelity loss 
can be separated into intrinsic errors, which include er- 
rors resulting from the use of weak-coupling and weak- 
driving approximations, and errors resulting from deco- 
herence [la, 113 ■ The latter, which of course depends 
sensitively on the experimental architecture and noise 
sources, is largely a function of total gate time. There- 
fore, by evaluating the intrinsic fidelity as a function of 
total gate time (which indirectly determines the strength 
of the qubit-qubit interaction), we can separate intrin- 
sic and extrinsic errors in a way that allows application 
to a wide variety of architectures and environments. The 
CNOT fidelity curves we calculate are the intrinsic fideli- 
ties as a function of gate time, optimized over all pulse 
sequence parameters and coupling constants that lead to 
the same total gate time. We calculate fidelity curves for 
coupled-qubit models with commonly occurring forms of 
interaction, such as XY and Heisenberg couplings, as well 
as for randomly generated ones with lower symmetry. We 
find that for qubit-qubit interactions not too close to that 
of the Ising model, the fidelity curves are largely insen- 



sitive to the form of the interaction. This allows us to 
provide a single fidelity curve for non-Ising-like models 
and a criterion for determining its applicability. 

The remainder of the paper is organized as follows: In 
Secim we review the perturbative CNOT gate design de- 
rived in Ref. [isj and describe the model considered there. 
In Sec. mil we discuss the various sources of intrinsic er- 
rors and define the fidelity measure used for all subse- 
quent computations. In Sec. IIV| we explain the method- 
ology used for our fidelity calculations, which involve 
exact numerical simulations of the underlying coupled- 
qubit models. In Sec.|Vl we present fidelity curves for dif- 
ferent forms of interaction and the interaction-averaged 
fidelity, and in Sec. IVII we explain the poor performance 
when the interaction is close to that of the Ising model. 



II. CNOT PROTOCOL 

In this section, we briefly review the main results of 
Ref. J^j for constructing a CNOT gate pulse sequence. 

A. Model Hamiltonian 

The Hamiltonian for a wide variety of physical sys- 
tems being considered for quantum computation can be 
written as 
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where J^j, is a 3 x 3 coupling matrix which takes differ- 
ent forms for different architectures under consideration. 
The parameters ci and f2i (with VLi ^ e^) are tunable and 
weak coupling {\Jfiu \ ^ q) is assumed. The CNOT gates 
are implemented according to a pulse sequence consisting 
of two entangling operations along with single qubit ro- 
tations. The entangling operations are carried out with 
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tuned qubits (ei = 62) and the local rotations are per- 
formed with detuned qubits. Table |T] gives values of the 
model parameters used for our simulations. For weakly 
coupled tuned qubits, the Hamiltonian(IT]) can be written 
in the interaction picture (or rotating frame) as 



(2) 



where 



The parameters J and J' in Eq.Q are given by 



(3) 



J = 



Jxx ~i~ '-^y 



and J = 



J, 



yx 



(4) 



TABLE I: Parameter values used in this work. The tuned 
qubit frequency is the frequency of the qubits used during 
the entangUng operations. During single-qubit operations, 
the frequency of qubit 2 is suddenly increased to 11 GHz. 
The ranges of allowed values of Rabi frequencies and overall 
coupling strengths are used to constrain the optimization de- 
scribed in the body of the paper. J*^ and g are defined in 
Eq. dm). 



parameter 


value 


common tuned qubit frequency e/h 
qubit-qubit detuning 
range of allowed Rabi frequencies Q,/h 
range of allowed coupling strengths g/h 

range of gate times tgatc considered 


10 GHz 
1 GHz 
50-500 MHz 
1-500 MHz 
10-50 ns 


Jp^ for Isotropic Heisenberg coupling 




"1 0" 

1 
1 




J^j, for Ising coupling 




"0 0" 


1 




J^i^^ for XY coupling 




"1 0" 
1 






B. CNOT pulse sequence 

The pulse sequence derived in Ref . fTsjl . carried out 
from right to left, is 



CNOT : 



Ry 



(5) 



where Rp,{0)i = e 2(^)'^i (with /i = x,y,z and i=l,2) is 
a single qubit rotation. Here 



If = arg (J + iJ') and At 



8V J2 + J'2 ■ 



(6) 



The operator ^--'^'HAt/h represents the action of bringing 
the qubits into resonance for a time At. The CNOT pulse 
sequence given in Eq. ([5]) involves two rotations about the 
z axis. For our exact simulations below, it will be conve- 
nient to rewrite ([5]) in terms of x and y rotations, leading 
to 



CNOT = e'- 



Ry{v^)2R^ (f). 



xe-^^^'/^R, {n), e-«^*/'' [r, (|) ^ R, (-|) 



(7) 

This is the CNOT pulse sequence used in the present 
analysis. Operations inside square brackets can be per- 
formed simultaneously. 



III. INTRINSIC FIDELITY 

The CNOT pulse sequence ([7|) is an exact identity; the 
errors come from realizing the individual terms in it using 
Hamiltonian ([T]). Here we discuss the possible sources of 
errors and the precise definition of fidelity used in this 
work. 



A. Sources of error 

As mentioned earlier, we are not concerned here with 
extrinsic errors, such as noise and decoherence, since 
these depend very much on the specific experimental ar- 
chitecture and noise sources. The fidelity loss computed 
here originates from intrinsic sources. The exact CNOT 
pulse sequence ^ is derived in the context of the approx- 
imate Hamiltonian which is derived from the model 
Hamiltonian ([IJ assuming weak coupling {\Jfj.u\ ^ Ci) 
and weak driving (fli <^ ei). These approximations con- 
tribute to the accumulation of fidelity loss. In addition, 
we also assume that the coupling is always on, even when 
the local rotations are being performed. This assump- 
tion is necessitated by the fixed coupling used in most 
experimental architectures. Due to the presence of more 
local rotations than entangling operations in the pulse 
sequence, the latter causes a larger contribution to the 
intrinsic error. 



B. Definition of fidelity 

In general, fidelity gives a measure of how close two 
quantum states are. There exists different measures of 
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fidelity. The definition we adopt is given by (See Ref.[l| 
Page 222, Eq.l4) 



F{\i^),p)^{^\p\ 



(8) 



where is considered to be a pure state and p is the 
density matrix of an arbitrary state. Here we are inter- 
ested in calculating fidelity between two operations (ideal 
CNOT and realized CNOT) which requires some modifi- 
cations of the definition given by Eq.®. In this context 
fidelity means how close these operations are. Suppose, 
we have two unitary operations U and ?7targct and we 
want to calculate the fidelity between these operations. A 
natural way to understand this closeness is to take a ran- 
domly generated vector \x) (defined on a Hilbert space), 
then apply the operations U and C/targot on the vector 
to obtain transformed vectors U \x) and C/taigct |x) : ^-^d 
finally identify these transformed states with p and |'0) 
in Eq.® to derive an expression for fidelity that depends 
on the state |x), 



Fx(C/targct,;7)=_(x|t/.^ 
1^1 



target (^IX>(x|f/^,.f/target|x); (9) 



Finally, we average over randomly generated |x) (cho- 
sen from a uniform distribution) to introduce an average 
fidelity, according to 

-^average (t/targct, U) = ]y"^j^ X/ (f^targot, U) , (10) 



where iV (|x)) is the total number of randomly generated 
Ix) states. To obtain a closed form expression of fidelity 
we change this sum to an integral, 



average 



(t^argct, U) 



MM\xWdv, 



(11) 



target ^ 

It has already been proven jl9l - l21| that, for any linear 
operator M on an n-dimensional complex Hilbert space. 



|(,|M|,)|^<iV-.?<™'' + l^<"'l 



52r, 



n{n + 1) 



(12) 



where the normalized state vectors \x) are defined on 
the unit sphere S'^""! in C. Using Eq.([12l) for a 4- 
dimensional Hilbert space, we can rewrite our expression 
for average fidelity as 



average (^target 



Tr(C/LgetC^) 



20 



(13) 



We use (ITS)) for computing fidelity between any two uni- 
tary quantum operations and express it in percent. 



IV. SIMULATIONS 

For a given qubit-qubit coupling tensor J^i,, the pulse 
sequence ([7]) is realized by performing the pair of two- 
qubit entangling operations with tuned qubits [for a time 
At given in ([6])] and the single-qubit operations with 
strongly detuned qubits. The time to implement the full 
pulse sequence depends on J^i, and the Rabi frequencies, 
which in principle can be different for each qubit and for 
each local rotation required. However, in this work we 
choose all Rabi frequencies to have the same value. 

The coupling tensor J^i, can be decomposed according 

to 



J^1V — 9 ^ '^uu: 



(14) 



where g > is a measure of the overall strength and 
J*y describes the form of the coupling. J*^ is defined to 
satisfy 



I J* J < 1 



all /i, ly. 



(15) 



Three important examples of J*^ are given in Table H] 

The simulations reported here are obtained by ex- 
act numerical integration of the model ([TJl.Our choice of 
fixed experimental parameters was motivated by super- 
conducting architectures [22| . We optimize over Rabi fre- 
quencies and qubit-qubit interaction strengths only, and 
do not allow for variation in the local rotation angles 
of Eq. ■ Although small refinements of these rotation 
angles can make slight improvements in the fidelity (by 
compensating for the qubit coupling that is suppressed 
by detuning but still always present), the fidelity changes 
are small on the scale of the main effects we consider 
(the dependence on the form of qubit-qubit interaction). 
These considerations lead us to vary the coupling tensor 
Jyjy, total gate time and Rabi frequency and compute fi- 
delity as a function of these variables. But since we know 
the pulse sequence, we can determine the total gate time 
as a function of J^^ and the single Rabi frequency by 
adding up the time required for each operation. So, for 
the simulation we fix total gate time, vary Rabi frequency 
within an allowed range given in Table U compute corre- 
sponding values of g and optimize the fidelity from the 
evolution of the original Hamiltonian((T|) . This procedure 
leads to a single point on a fidelity curve. 



V. FIDELITY CURVES 

Figures [T] and [5] give the optimal CNOT fidelities as 
a function of total gate time for the Heisenberg and XY 
interactions and Tables IhI and IIIII show corresponding op- 
timal values of relevant parameters. The fidelity curves 
are similar, indicating that a fidelity of 99% can be ob- 
tained in less than 15 ns and 99.9% can be obtained in 
about 50 ns. Alternatively, these results indicate that 
for these common forms of qubit-qubit coupling, 99% 
can be achieved with a coherence time in excess of 15 
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ns and 99.9% can be achieved with at least 50 ns of co- 
herence. [We remind the reader that the model ([H does 
not include higher energy (non-qubit) states, which fur- 
ther limit performance, and that results are obtained for 
10 GHz qubits.] 



TABLE II: Optimum values of parameters for Heisenberg in- 
teraction. 







n/h ('MU7\ 

\j 1 1 ii y ivj.±±zj I 


Qlh IMV{7\ 

i L 1 1 L I 1VJ.J.±Zj j 


10.00 


97.8321 


19.1964 


430 


11.25 


98.4599 


16.1049 


430 


12.50 


98.8405 


13.8710 


430 


13.75 


99.0881 


12.1813 


430 


15.00 


99.2579 


10.8586 


430 


16.25 


99.3792 


9.7950 


430 


17.50 


99.4688 


8.9212 


430 


18.75 


99.5368 


8.1905 


430 


20.00 


99.5895 


7.5704 


430 


22.50 


99.6646 


6.5749 


430 


25.00 


99.7144 


5.8108 


430 


27.50 


99.7489 


5.2058 


430 


30.00 


99.7794 


4.8851 


340 


40.00 


99.8452 


3.5124 


340 


50.00 


99.8734 


2.7419 


340 



interactions, rj = 1. A "typical" value of rj, defined by 
averaging the function 'ri{J^^) over an unconstrained uni- 
form distribution of J^u tensors, is about 0.52. 




97.5 



10 20 30 40 

total gate time (in nanoseconds) 



50 



FIG. 1: (Color online) Plot of optimal fidelity versus total 
CNOT gate time for the Heisenberg interaction. 
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TABLE III: Optimum values of parameters for XY interac- 
tion. 



total time (ns) 


fidelity [%] 


g/h (MHz) 


Q./h (MHz) 


10.00 


98.1750 


17.8571 


500 


11.25 


98.8618 


23.8095 


250 


12.50 


99.2710 


19.2308 


250 


13.75 


99.4902 


16.6667 


240 


15.00 


99.6174 


14.2857 


240 


16.25 


99.6966 


12.5000 


240 


17.50 


99.7494 


11.1111 


240 


18.75 


99.7864 


10.0000 


240 


20.00 


99.8133 


9.0909 


240 


22.50 


99.8491 


7.6923 


240 


25.00 


99.8713 


6.6667 


240 


27.50 


99.8861 


5.8824 


240 


30.00 


99.8973 


5.2083 


250 


40.00 


99.9211 


3.6765 


250 


50.00 


99.9311 


2.8409 


250 



The weak dependence of the fidelity curve on the form 
of interaction is typical, unless the interaction is close to 
that of the Ising model (see Table H]). To quantify this 
closeness we define a parameter [recall ^ and ([T4| ] 



■q \ 



V J2 + J'^ 
9 



(16) 



It can be shown that < 77 < \/2. For the Ising in- 
teraction, 77 = 0, whereas for the Heisenberg and XY 
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FIG. 2: (Color online) Optimal fidelity versus total gate time 
for the XY interaction. 

In Fig. 13] we show fidelity curves for randomly gener- 
ated forms of interaction with three fixed values 77. The 
unambiguous loss of fidelity for 77 = 0.1 is due to the fact 
that the interaction is close to the Ising limit (77 = 0). 
The similar behavior of fidelity for rj — 0.5 and 77 = 1.0 
affirms our assertion that fidelity curves do not signifi- 
cantly depend on the form of the interaction as long as 
77 is not too close to zero. The reason for the poor per- 
formance for small 77 is discussed in Sec. IVIl 

Given that the fidelity is largely independent of the 
form of interaction, as long as rj is not too small, it 
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10 20 30 40 50 



total gate time (in nanoseconds) 

FIG. 3: (Color online) Plot of optimal fidelity versus total 
gate time for random interactions. 



is useful to average over interaction forms to obtain 
interaction-independent fidelity curves. This is provided 
in Fig. m wliicli present interaction-averaged fidelity 
curves for 77 = 0.1,0.5, and 1.0. 




10 20 30 40 50 



total gate time (in nanoseconds) 

FIG. 4: (Color online) Plot of optimal fidelity versus to- 
tal gate time, averaged over randomly generated interactions 
with fixed values of rj, for = 0.1, 0.5, and 1.0. 



VI. CONCLUSIONS 

We have shown that the intrinsic fidelity versus total 
gate time for the CNOT gate as constructed by Eq. ([7]) 
is largely independent of the form of qubit-qubit interac- 
tion, as long as that interaction is not too close to that 
of the Ising model, as measured by the the parameter 



77 defined in (|16p . For typical (non-Ising-like) couplings, 
the fidelity is given in Fig. here one can use either the 
rj = 0.5 or i] — 1.0 curves. 

TABLE IV: Optimum values of parameters for Heisenberg 
interaction with 500 ns. amplitude damping. 



total time (ns) 


fidelity [%] 


n/h (MHz) 


fl/h (MHz) 


10.00 


96.9272 


21.0227 


370 


11.25 


97.3098 


17.3709 


370 


12.50 


97.4493 


14.8000 


370 


13.75 


97.4699 


12.8920 


370 


15.00 


97.4260 


11.4198 


370 


16.25 


97.3423 


10.2493 


370 


17.50 


97.2355 


9.2965 


370 


18.75 


97.1125 


8.5057 


370 


20.00 


96.9781 


7.8390 


370 


22.50 


96.6904 


6.7766 


370 


25.00 


96.3866 


5.9677 


370 


27.50 


96.0740 


5.3314 


370 


30.00 


95.7565 


4.8177 


370 


40.00 


94.4736 


3.4774 


370 


50.00 


93.1969 


2.7206 


370 




10 20 30 40 50 

total gate time (in nanoseconds) 



FIG. 5: (Color online) Plot of optimal fidelity versus total 
gate time for Heisenberg interaction in presence of amplitude 
damping with 500 ns. 

The origin of the poor fidelity when rj is small can be 
understood as follows: In the pulse construction ([7]) of 
Ref. [l3|, a Cartan decomposition is used to decompose 
the time-evolution operator generated by ^ into single- 
qubit rotations, an entangling operator, and a global 
phase factor. The entangler has the form 

A{x,y,z) EE e'i{'-^<^f<^2+y<^1<^l+-<-f<^S)^ (17) 

where x, y, and z are three coordinates (angles). Fol- 
lowing Zhang et al. Q, the entangler coordinates trace 
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out a trajectory in the three-dimensional space of entan- 
glers as time progresses. In the construction of Ref. [l3| . 
the trajectory starts in the plane x = y, and then a re- 
focusing TT pulse is used to reflect the trajectory to the 
point (|-,0, 0) or (— 1,0,0) on x axis. (The actual point 
reached depends on the sign of J.) The time it takes to 
do this — neglecting the time needed for the tt pulse — is 
2 At [see or irh/Agrj. Including all the single-qubit 
rotations in ([7]) leads to a total gate time of 



Agri 



n 



(18) 



Because the first term in (|18l) is inversely proportional 
to Tjg, for a fixed gate time a larger value of coupling 
strength g is required when r] is small. But when 
increases the assumption of weak coupling used in \lc 
is violated and the corrections to the rotating-wave- 
approximation get larger. Furthermore, that large cou- 
pling leads to considerable errors during the single-qubit 
operations because the qubit-qubit interaction is not 
switched off. 

Although we have focused on the intrinsic fidelity in 
this work, it is interesting to calculate one example of a fi- 
delity curve with decoherence. We choose the Heisenberg 
interaction for this study, with 500 ns amplitude damp- 



ing. Reoptimizing n and g for each total gate time leads 
to the fidelity curve shown in Fig.®. Table ITVl gives the 
corresponding optimal parameters. The curve exhibits a 
maximum fidelity (« 97.47%) at about 13.75 ns. which 
represents the optimal time to construct a CNOT with 
this assumed decoherence model. The optimal interac- 
tion strength and Rabi frequency are g/h w 12.89 MHz 
and n/h^ 370 MHz. 

Figs. ([3]) and ^ are our principal results. To use the 
fidelity curve of Fig.Q for a particular application, one 
should calculate the 77 value for the application and ex- 
trapolate between the curves provided. We note, how- 
ever, that for small 77 the pulse sequence ([7]) is not useful, 
and one should construct an alternative pulse sequence 
using the methods of Refs. Q and jll] to generate an 
entangler on the z axis instead of the x axis. 
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